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Abstract 

Based on the vector angular spectrum method and the stationary phase 
method and the fact that a circular aperture function can be expanded into 
a finite sum of complex Gaussian functions, the analytical vectorial structure 
of a four-petal Gaussian beam (FPGB) diffracted by a circular aperture is de- 
rived in the far field. The energy flux distributions and the diffraction effect 
introduced by the aperture are studied and illustrated graphically. Moreover, 
the influence of the /-parameter and the truncation parameter on the non- 
paraxiality is demonstrated in detail. In addition, the analytical formulas 
obtained in this paper can degenerate into un-apertured case when the trun- 
cation parameter tends to infinity. This work is beneficial to strengthen the 
understanding of vectorial properties of the FPGB diffracted by a circular 
aperture. 
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1. Introduction 

In recent years, there have been increasing interests in the study of beam 
pattern formation. Many beam patterns have potential and practical ap- 
plications. For example, many researches have been done on dark-hollow 
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beam because it is a powerful tool in precise manipulation and control of 
microscopic particles [II, 0, 0, 0, 111 . In fact, many other kinds of special laser 
patterns, such as flower-like patterns and daisy patterns have been observed 
and investigated js), 0|. More recently, a new form of laser beam called 
four-petal Gaussian beam has been introduced and its properties of pass- 
ing through a paraxial ABCD optical system have been studied M. Subse- 
quently, much work has done on four-petal Gaussian beam 0,0,0,0,113 



As is well known, some researches and applications are conducted in the far 
field. Meanwhile, in practical applications, however, there are more or less 
aperture effects, so it is of practical significance to study the influence of a 
hard-edged aperture on the far-field properties of four-petal Gaussian beam. 

In this paper, firstly the power transmissivity of the truncated four-petal 
Gaussian beam passing through a circular aperture is studied. Secondly, the 
analytical vectorial structure of four-petal Gaussian beam diffracted by a 



circular aperture is derived based on vector angular spectrum method [14 



• ap 

15l . Il6l Il7| . stationary phase method 18], and the fact that a circular aper- 



ture function can be expanded into a finite sum of complex Gaussian func- 
tions [l9[. Based on the analytical vectorial structure of four-petal Gaussian 
beam diffracted by a circular aperture, the energy flux expressions of TE 
term, TM term and the whole beam are also obtained, respectively. Thirdly, 
some typical numerical examples are given to illustrate the influence of the 
diffraction effect introduced by an aperture on the far-field energy flux dis- 
tributions of four-petal Gaussian beam. Furthermore, the influence of the 
/-parameter and the truncation parameter on the non-paraxiality is demon- 
strated. Finally, some simple conclusions are given. 



2. Analytical vectorial structure of a hard-edged-diffracted four- 
petal Gaussian beam 

Let us consider a half space z > filled with a linear, homogeneous, 
isotropic, and nonconducting medium characterized by electric permittivity 
e and magnetic permeability [i. All the sources only lie in the domain z < 0. 
For convenience of discussion, we consider a four-petal Gaussian beam with 
polarization in the x direction, which propagates toward the half space z > 
along the z axis. The transverse electric field distribution of the incident 
four-petal Gaussian beam at z = plane can be written by [s| 
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E x (x,y,0) = G n (^} e xp(-^-^,n = 0,l,2..., (1) 

E v (x,y,0) = Q, (2) 

where G n is a amplitude constant associated with the order n; u>o is the 
1/e 2 intensity waist radius of the Gaussian term; n is the beam order of the 
four-petal Gaussian beam. When n — 0, Eq. p]) reduces to the ordinary 
fundamental Gaussian beam with the waist being wq at the plane z = 0. 
The time factor exp(— iuit) has been omitted in the field expression. After 
simple calculation, we know that the distance of diagonal petals is given by 

d = 2(2n) 1/2 w . (3) 

According to Eq. [31 the distance d is determined by beam order n and waist 
size wo- 

Supposing that a hard-edged circular aperture is located at the z = 
plane, and the center of the circular aperture is the origin. The geometry 
for the screen with a circular aperture and the coordinate system is shown 
in Fig. [TJ The corresponding circ function can be written as follows 

circ(^) = { P ~ (A) 
V \ 0, p>R. [V 

where R denotes the radius of the circular aperture and p = \J x 2 + y 2 is ra- 
dial distance. The power transmissivity of the truncated four-petal Gaussian 
beam is given by 



T 



J o R C\E x (p,e,0) \ 2 pdpd6 



n roo p2tt 



f™f 2% \E x (p,e,0) \ 2 pdpd9 
T(4n + l,2/3 2 ) 
4nr(4n) ' 



(5) 



where T(z) = J Q tf^e^dt is Euler gamma function, and T(a, z) = f z t a ~ x e~ l dt 
is incomplete gamma function. = R/wq is defined as the truncation param- 
eter. When n tends to zero, one can also obtain the fundamental Gaussian 
beam power transmissivity, namely 

T = l-exp(-2i?7^o)- (6) 
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Figure 1: Illustrating the geometry of the screen with a circular aperture and the coordi- 
nate system. 
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Figure 2: The power transmissivity of the truncated FPGB versus the truncation param- 
eter /3 for various beam order n = (dash-dotted curve), n = 1 (solid curve), and n = 3 
(dashed curve), respectively. 



Fig. [2]represents the power transmissivity of the truncated FPGB versus the 
truncation parameter for various beam order n = (dash-dotted curve), 
n = 1 (solid curve), and n = 3 (dashed curve), respectively. It is found that 
T n increases rapidly with decreasing n for the same value of (3. 

Starting from Maxwell's equations, the basic principles of vector angu- 
lar spectrum method are vectorial plane wave expansion. The propagating 
electric field toward half space z > turns to be [14|, [l7l |20 



r poo 

E x (x,y,z)= / A x (p,q,7)exp[ik(px + qy + ^z)]dpdq, (7) 

J J —OO 

p poo 

E y (x,y,z)= / A y (p,q,"/)exp[ik(px + qy + jz))dpdq, (8) 



E z (x,y,z) 




P Q 
-A x (p,q,j) + -Ay(p,q,"/) 

7 7 
x exp[ik(px + qy + •yz)]dpdq, 



(9) 



where k = 2n/\ denotes the wave number in the medium related wave length 
A, 7 = a/1 — p 2 — q 2 , if p 2 + q 2 < 1 or 7 = %\/p 2 + q 2 — 1, if p 2 + q 2 > 1. The 
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values of p 2 +q 2 < 1 correspond to the homogeneous waves which propagate at 
angles arcsin ^/p 2 + q 2 with respect to the positive z axis, whereas the values 
of p 2 + q 2 > 1 correspond to the evanescent waves which propagate along 
the boundary plane but decays exponentially along the positive z direction. 
In terms of Fourier transform, the transverse components of the vectorial 
angular spectrum of the electric field just behind the aperture are expressed 
as follows 



Iff 00 p 
A x(p,q) = -^J J E x (x,y,0)ciTc(—)exp[-ik(px + qy)\dxdy, (10) 



A y{p,q) 



1 

Y 2 



r r°° p 

/ / Ey(x, y, 0)circ(— ) exp[— ikypx + qy)]dxdy. (11) 

J J-oo R 



As is well known, circ function can be expanded into a finite sum of complex 



Gaussian functions |19| . |2lL |22 



N 



chc(— ) = ^Aiexp(-Bi 



i=i 



R 2 ' 



;i2) 



where Ai, Bi are the coefficients and iV is the number of complex Gaussian 
terms; they can be found in Table 1 of Ref 19(. On substituting Eqs. ([T]), 
© and (H2D in Eqs. ([TP]) and (fTTll. one can find 



A x (p,q) 



G, 



4vr 2 / 2 



r(n + -) 



2 N 



1 1 

x i F i I n + o; ~; 



J2 A 

h 
P 



i=i 



P 2 + B l 



2n+l 



2' 2' 4/ 2 /3 2 + 5; 



1 1 

Xi F i I ™+ tt; 



q' 



0' 



(13) 



2' 2' Af 2 /3 2 + B t 

A y (p,q) = 0, (14) 

where/ = 1/kwo, which is /-parameter, and i-Fi(-) denotes confluent hy- 
pergeometric function. It is well known that Maxwell's equations can be 
separated into transverse and longitudinal field equations and an arbitrary 
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polarized electromagnetic beam, which is expressed in terms of vector an- 
gular spectrum, is composed of the transverse electric (TE) term and the 
transverse magnetic (TM)term, namely, 

E(r) = E TE (r)+E TM (r), (15) 



where 



H(r) = H TE {r) + H TM (r), (16) 



1 

p2 _|_ g2 

x exp(iku)dpdq, (17) 




E T E(r) =11 2 ~ 2 [qA x (p, q) - pA y (p, q)] (qe x - pe y ) 



OO 

H T E(r) = J- I I 2 2 [qA x (p, q) - pA y (p, q)]{jrfe x + qye y - b 2 e z ) 
V A 4 J J -oo P ~r q 

x exp(iku)dpdq, (18) 



and 



-> f f°° 1 b 2 
E T m{t) = 2 j 2 [p4r(p, g) + 54/ (p, 9)](pe x + ge y e z ) 

x exp(iku)dpdq, (19) 



H T M{r) = ~]J~ l J J [p4(p»9) + q A y(p,q)}-^(qex-pe y ) 

x exp(iku)dpdq, (20) 

where r = xe x +ye y +ze z is the displacement vector and e x , e y , e z denote unit 
vectors in the x, y, z directions, respectively; u = px + qy + •yz; b 2 = p 2 + q 2 . 

Generally speaking, the evolution of beam is often studied by virtue of 
numerical simulation. However, in the far field framework, the condition 
k(x 2 + y 2 + z 2 ) 1 / 2 — > oo is satisfied due to z is big enough. By virtue of 
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the method of stationary phase 18j|, the TE mode and the TM mode of the 
electromagnetic field can be given by 



E TE (r) 



.G n Z R yz 

' 2~~2 

7r r z pr 
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Xl 1 V n+ 2 ; 2 ; 4/ 2 r 2 /3 2 + Z^ 
x exp(ikr)(ye x - xe y ), 



(21) 



^ G n Z R yz 
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(22) 



and 



E T M(r) 



■ G n Z R x 
7r r z p z 
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B 2 
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1 1 
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y 2 /? 2 



2' 2' 4/ 2 r 2 /3 2 + 5, 
x exp(iA;r)(xze x + yze y — p 2 e z ), 



(23) 
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2' 2' AfWpt + Bi 
x exp(ikr)(ye x - xe y ), 



(24) 



where r = \/ x 2 + y 2 + z 2 , and Zr = ixw 2 /\ is Rayleigh length. Eqs. (J2il — (|24l) 
are analytical vectorial expressions for the TE and TM terms in the far field 
and constitute the basic results in this paper. It follows that spherical wave 
front remain unchanged for apertured FPGB in the far field. The results 
obtained here are applicable for both non-paraxial case and paraxial case. 
From Eqs. (1211) — ff24|) . one can find that 

E TE {r) ■ E TM (r) = 0, (25) 

H TE (r) ■ H TM {r) = 0. (26) 

According to Eqs. ( |25i) and ( |26i) . the TE and TM terms of a hard-edged- 
diffracted four-petal Gaussian beam are orthogonal to each other in the far 
field. 

3. Energy flux distributions in the far field 

The energy flux distributions of the TE and TM terms at the z = const 
plane are expressed in terms of the z component of their time-average Poynt- 
ing vector as 

(S z )te = \Re{E TE x H TE ] Z , (27) 

(S z )tm = -^R e \E T M x Htm}z, (28) 

where the Re denotes real part, and the asterisk denotes complex conjuga- 
tion. The whole energy flux distribution of the beam is the sum of the energy 
flux of the TE and TM terms, namely 

(S z ) = (S Z ) TE + (S Z ) TM . (29) 
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On substituting Eqs. (12T|)- (J23j) in Eqs. (127j)- (125]) yields 



2 V /x 7r 2 r 5 p 2 



1 \eG\Z\ x 2 z 



-\ i 



2 V /i 7r 2 r 3 p 2 



r(n + -) 



r(n + -) 



S n (x,y,z, f,p) 



S n (x,y,z,f,(3) 



(30) 
(31) 



Therefore, the whole energy flux distribution of a hard-edged-diffracted four- 
petal Gaussian beam in the far field is given by 

1 fTG 2 Z 2 



where 



2 V [i it 2 r 3 p 2 
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the function 5 n (-) is defined as above for simplifying expressions of energy 
flux. Eqs. (1301)— ( 1321) indicate that the diffraction effect introduced by a cir- 
cular aperture is described by the truncation parameter j3 in the far field. The 
smaller the truncation parameter is, the more strongly the field is diffracted 
by the aperture. In addition, as the truncation parameter tends to infinity, 
Eqs. (SOD — d32D degenerate into 
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(35) 



10 



(Sm) 



1 ^ G 2 n Z\ z 

2 Y jj, 7r 2 r 3 p 2 




(36) 



As a matter of fact, Eqs. (134?) — fl36|) are energy flux expressions in un-apertured 
case, which is not difficult to understand. 

4. Numerical examples and discussion 

The normalized energy flux distributions of the TE term, the TM term 
and the whole of the un-apertured beam at the plane z = 2000A for different 
beam order n = 1 and 3 versus x/X and y/X are illustrated in Fig. [3] based on 
Eqs. ( 134)) — ( 1361) . respectively. Waist size wo is set to 10A in the calculation. 
Apparently, the TE term and the TM term are orthogonal to each other. 
The four-petal Gaussian beam splits into a number of small petals in the 
far field, which differs from its initial four-petal shape. Furthermore, the 
number of petals in the far field gradually increases when the parameter n 
increases, which has potential applications in micro-optics and beam splitting 
techniques, etc jsj. The above conclusion is also applicable to four-petal 
Gaussian beam diffracted by a circular aperture. In fact, the four-petal 
Gaussian beam with beam order n is not a pure mode, which can be regarded 
as a superposition of n 2 two dimensional Hermite-Gaussian modes |8|, and 
different modes evolve differently within the same propagation distance. The 
overlap and interference between different modes result in the propagation 
properties of the four-petal Gaussian beam in the far field. 

The influence of the diffraction effect introduced by the aperture on the 
energy flux distributions of the four-petal Gaussian beam is depicted in Fig. @] 
and Fig. based on Eqs. (15U1) — (1521 . Waist size w is set to 10A. For 
simplicity, the beam order n is set to be 1. In Fig.lH the truncation parameter 
is set to ft — a/2 and (3 = 1, respectively. According to Eq.(j5J, the four peak- 
value positions of the incident beam is just on the boundary line of the circle 
when the truncation parameter is set to ft = a/2 for n = 1. Comparing 
top subfigure of Fig. [5J with Fig. HI one can find that central spot and side 
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Figure 3: (Color online) Normalized energy fluxes (S z )te, (S z )tm and (S z ) of the FPGB 
with beam order n = 1 and n = 3 (from top to bottom) in the reference plane z = 2000 A, 
respectively. Waist size wq is set to 10A, and the circular aperture radius R — > oo. 
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Figure 4: (Color online) Normalized energy fluxes (S z )te, (S z )tm and (S z ) of the hard- 
edged-diffracted four-petal Gaussian beam with beam order n = 1 in the reference plane 
z = 2000 A. The truncation parameter is set to /3 = y/2 and (3 = l(from top to bottom), 
respectively. Waist size wo is set to 10A. 
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Figure 5: (Color online) Cross section and on-axis value of the normalized total energy 
flux of the hard-edged-diffracted four-petal Gaussian beam with the beam order n = 1 in 
the reference plane z = 2000A. Waist size wo is set to 10A. (a) The cross section with 
respect to x direction (y = 0) with various truncation parameter ft = 1 (dashed curve), 
ft = \/2 (circled curve), and ft — oo (solid curve), respectively, (b) Normalized axial energy 
flux at z = 2000A versus the truncation parameter ft. 

lobes spread more widely in the far field when the circular aperture exists. 
Moreover, the smaller radius of the aperture is, the more wide distribution 
of the energy flux is. This phenomenon is easy to understand because the 
initial field is confined by the aperture and its diffracted field has bigger 
divergence angle in the far field [23|. For the sake of showing clearly, cross 
section and on-axis value of the whole energy flux is plotted based on Eq. ( 1321) 
in Fig. El Fig. EJa) shows cross section of the total energy flux at y = 
for different truncation parameter, from which we can see clearly that the 
relative values of the lobes and the full width at half maximum (FWHM) of 
central spot become bigger when the truncation parameter decreases. The 
values of FWHM are 41.6A for R = oo, 61.8A for = ^2, and 83.8A for 
(3 = 1, respectively. In addition, the ratios of the maximum value of the first 
side-lobe to that of the central spot are 0.19 for /3 = oo, 0.32 for /3 = \/2, 
and 0.34 for (3 = 1, respectively. The on-axis energy flux of the whole beam 
at z = 2000A versus truncation parameter /3 is shown in Fig. 15(b). It should 
be noticed that truncation parameter cannot be too small in order to obtain 
the better transmissivity in practical application. 

It is well known that the paraxial approximation is allowable when wq/ A 
is much larger than 1. The far field expression of the four-petal Gaussian 
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Figure 6: (Color online) Cross section at y = of the normalized total energy flux of 
the hard-edged-diffracted four-petal Gaussian beam with the beam order n = 1 in the 
reference plane z = 2000 A. The solid and the dashed curves denote the non-paraxial 
and the paraxial results, respectively. Truncation parameter /3 is 2. In subfigure(c) cross 
section at x — is also plotted based on non-paraxial result, (a) / = 0.1, (b) / = 0.2, (c) 
/=1.5. 
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Figure 7: (Color online) Cross section at y = of the normalized total energy flux of 
the hard-edged-diffracted four-petal Gaussian beam with the beam order n = 1 in the 
reference plane z — 2000A. The solid and the dashed curves denote the non-paraxial and 
the paraxial results, respectively. The f-parameter is set to be 0.1. (a) j3 = 1.1, (b) (3 — 1.5, 
(c) /3 = 2.5. 
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Figure 8: (Color online) Normalized total energy flux of the hard-edged-diffracted four- 
petal Gaussian beam with the beam order n = 1 in the reference plane z = 2000A. The 
/-parameter is set to be 1.5. /? = 2.0. (a) non-paraxial result, (b) paraxial result. 



beam diffracted by a circular aperture under paraxial regime can be treated 
as a special case by using approximation 

P 2 

r = z + y ~ z - ( 37 ) 

Therefore, the paraxial expression of the whole energy flux of the four-petal 
Gaussian beam diffracted by a circular aperture turns out to be 

4 

\S n (x,y,z,f,P) | 2 . (38) 

Eq. (1381) is symmetric about the x and y variables, whereas Eq. ( 1321) is some- 
what asymmetric about the x and y variables. With the given beam order, 
the non-paraxiality of an apertured FPGB depends on /-parameter and trun- 
cation parameter /3. Fig. [6] shows cross section at y — of the normalized 
total energy flux of the hard-edged-diffracted four-petal Gaussian beam with 
different /-parameter in the reference plane z = 2000 A. Truncation param- 
eter P is 2. For simplicity, the beam order n = 1 (hereafter). The solid and 
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the dashed curves denote the non-paraxial and the paraxial results, respec- 
tively(hereafter). From Fig. |?^a) one find that the difference between the 
non-paraxial and the paraxial results is negligible for /3 = 2 and / < 1. 
The difference between them becomes evident as /-parameter increases. The 
central beam spot obtained by paraxial result is obviously larger than that 
obtained by non-paraxial result and side lobes disappear when /-parameter 
increases enough, which is shown in Fig. El^c). Moreover, the cross section 
at x = is also plotted in Fig. [6]^c) in order to observe the asymmetry. 
Fig. [7] shows cross section at y = of the normalized total energy flux of 
the hard-edged-diffracted four-petal Gaussian beam with different truncation 
parameter in the reference plane z = 2000A. The /-parameter is 0.1. The 
difference between paraxial result and non-paraxial result becomes evident as 
truncation parameter decreases. Comparing Fig. [6] with Fig. we obtained 
a conclusion that the /-parameter plays a more key role in determining the 
non-paraxiality of an apertured FPGB than does the truncation parameter 



(3, which is similar to the conclusions of Ref. [24], [25(. In order to under- 
stand, Fig. E^c) are plotted in contour graphs, which is shown in Fig. [SJ It 
can be clearly shown the non-paraxial result is approximately elliptic, which 
is obviously different from the paraxial result. The above conclusion is also 
applicable to other higher beam order. 



5. Conclusions 

In summary, the vectorial structure of an apertured four-petal Gaussian 
beam in the far field is derived in the analytical form by using the vector 
angular spectrum method, the complex Gaussian expansion of the circular 
aperture function, and the stationary phase method. Based on the analyti- 
cal vectorial structure of an apertured beam, the energy flux distributions of 
the TE term, the TM term and the whole beam are derived in the far-field. 
Our formulas obtained in this paper are applicable to both non-paraxial case 
and paraxial case. When the truncation parameter tends to infinity, our 
formulas degenerate into the un-apertured case. The four-petal Gaussian 
beam cannot preserve its initial shape, and the number of petals in the far 
field gradually increases when beam order n increases. Energy distributions 
spread more widely in the far field when the circular aperture exists. The 
influence of /-parameter and truncation parameter on the far-field behavior 
is also studied in detail. The /-parameter plays a more key role in deter- 
mining the non-paraxiality of an apertured FPGB than does the truncation 
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parameter /3. In addition, the asymmetry of beam spot becomes apparently 
with increasing non-paraxiality. This work is important to understand the 
theoretical aspects of vector FPGB and is beneficial to its practical applic- 
tion. 
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